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ABSTRACT 

We study the cross-correlation between galaxies of different luminosities and colors, using a sample 
of 284489 galaxies selected from the Sloan Digital Sky Survey, Data Release 4. Galaxies are divided 
into 6 samples according to luminosity, and each of these samples is further divided into red and 
blue subsamples according to their (g — r) colors. Projected auto-correlation is estimated for each 
subsample, and projected cross-correlation is estimated for each pair of subsamples. At projected 
separations r p > 1 /i Mpc, all correlation functions are roughly parallel, although the correlation 
amplitude depends systematically on both luminosity and color. On r p < 1 ft. Mpc, the auto- and 
cross-correlation functions of red galaxies are significantly enhanced relative to the corresponding 
power laws obtained on larger scales. Such enhancement is absent for blue galaxies and in the cross- 
correlation between red and blue galaxies. We estimate the relative bias factor on scales r > 1 ft,~ 1 Mpc 
for each subsample using its auto-correlation function as well as its cross-correlation functions with 
other subsamples. The relative bias factors obtained from different methods are similar. For blue 
galaxies the luminosity-dependence of the relative bias is strong over the entire luminosity range probed 
(—23.0 < M r < —18.0), while for red galaxies the dependence is weaker and becomes insignificant 
for luminosities below L*. In order to examine whether a significant stochastic/nonlinear component 
exists in the bias relation, we study the ratio TZij = WuWjj/Wf^ where is the projected correlation 
between subsample H' and t j'. We find that the values of TZij ar e all consistent with 1 for all-all, red- 
red and blue-blue samples, however significantly larger than 1 for red-blue samples. For faint red - 
faint blue samples the values of TZij are as high as ~ 3 on small scales (r p < 1 /i _1 Mpc) and decrease 
with increasing r p . These results suggest that a significant stochastic/nonlinear component exists in 
the relationship between red and blue galaxies, particularly on small scales. 
Subject headings: large-scale structure of the universe - galaxies: halos - methods: statistical 



1. INTRODUCTION 

In the current scenario of structure formation, galaxies 
are assumed to form in the cosmic density field through 
a series of physical processes, such as non-linear grav- 
itational collapse, gas cooling and star formation. In 
this scenario the distribution of galaxies in space is ex- 
pected to trace the underlying density field to some de- 
gree, but the relationship is not expected to be perfect, 
because many of the processes involved in galaxy for- 
mation can complicate the relationship between galaxies 
and the dark matter density field. Thus, the study of 
galaxy clustering in space is an important step both in 
understanding the mass distribution in the universe, and 
in understanding how galaxies form in the cosmic density 
field. The statistical tool commonly adopted to quantify 
galaxy clustering in space is the correlation functions of 
galaxies (Peebles 1980). During the last two decades, 
many authors have estimated the two-point correlation 
functions of galaxies using various redshift catalogs, to 
study how galaxy clustering depends on the properties 
of galaxies, such as luminosity (e.g. Borner, Mo & Zhou 
1989; Park et al. 1994; Loveday et al. 1995; Guzzo et al. 
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1997; Benoist et al. 1996; Norberg et al. 2001; Zehavi et 
al. 2002; Zehavi et al. 2005; Li et al. 2006), color (e.g. 
Willmer et al. 1998; Brown et al. 2000; Zehavi et al. 
2002; Zehavi et al. 2005; Li et al. 2006), spectral type 
(Norberg et al. 2002; Budavari et al. 2003; Madgwick et 
al. 2003), morphological type (e.g. Jing, Mo & Borner 
1991; Guzzo et al. 1997; Willmer et al. 1998; Zehavi 
et al. 2002; Goto et al. 2003), stellar mass, stellar sur- 
face mass density and concentration (Li et al. 2006). All 
these analyses show that galaxies of different properties 
may have different distributions in space, indicating that 
galaxies are biased tracers of the underlying density field 
and that one has to be careful when using the observed 
galaxy distribution to infer the mass distribution in the 
universe. 

The observed correlation functions of galaxies can pro- 
vide important constraint on the relationship between 
galaxies and dark matter halos, as is demonstrated 
clearly in the recent halo-occupation and conditional- 
luminosity- function models (Jing, Mo & Borner 1998; 
Seljak 2000; Peacock & Smith 2000; Berlind & Weinberg 
2002; Cooray k, Sheth 2002; Yang, Mo & van den Bosch 
2003; van den Bosch, Yang & Mo 2003). With accurate 
measurements of the clustering properties for galaxies of 
different properties, it is now possible to study in great 
detail how different galaxies occupy different halos and 
how galaxies trace the cosmic density field. 

The correlation analyses carried out so far are mostly 
based on the auto-corrclation function of galaxies. Al- 
though such analyses can provide important informa- 
tion about how a population of galaxies is distributed 
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in space, they do not tell us anything about the rela- 
tionship between galaxies of different properties. For ex- 
ample, two populations of galaxies may each be strongly 
clustered in space and yet be spatially segregated so that 
the cross-correlation between them is weak. In reality, 
all populations of galaxies may be biased tracers of the 
underlying mass distribution, and so they must all be 
positively correlated in space. However, the amplitude 
of the cross-correlation between any two populations of 
galaxies depends not only on the bias factors of these 
two populations, but also on how well each of them trace 
the mass density field, i.e. on the stochasticity in the 
bias relation. Thus, it is important to measure the cross- 
correlation functions between galaxies of different prop- 
erties, so as to quantify the relationships between the 
spatial distributions of different galaxies. 

In this paper, we use galaxy samples constructed from 
the Sloan Digital Sky Survey Data Release 4 (SDSS- 
DR4) to measure both the auto-correlation functions for 
galaxies of different luminosities and colors and the cross- 
correlation functions between different galaxies. We com- 
pare the shapes of the various correlation functions to 
check the validity of the assumption of linear bias. For 
each population of galaxies (i.e. for galaxies within given 
ranges of luminosity and color), we estimate its rela- 
tive bias factors using its auto-correlation function as 
well all its cross-correlation functions with other popu- 
lations of galaxies. These bias factors are averaged to 
give a mean bias factor for each population, and com- 
pared among themselves to constrain possible stochas- 
tic/nonlinear component in the bias relation (Mo & 
White 1996; Dekel & Lahav 1999). This paper is or- 
ganized as follows. In section [21 we describe the galaxy 
samples to be used and how random samples are con- 
structed. In section [31 we describe our method for es- 
timating the correlation functions. Our results about 
the luminosity- and color-dependence of the correlation 
functions are presented in Sections 0] and respectively. 
In Section [6] we discuss the implications of our results 
for the possible existence of stochastic/no n- linear bias in 
galaxy distribution. Finally, we summarize our findings 
in Section [71 

Throughout this paper, galaxy distances are obtained 
from rcdshifts assuming a cosmology with Q m .o = 
0.3, f^A.o — 0.7. Distances are quoted in units of 
/i _1 Mpc, and absolute magnitudes are quoted in terms 
of M — 5\ogh, i.e. with h set to be 1, where h — 
ffoAlOOkms^Mpc- 1 ). 

2. THE DATA 

2.1. The SDSS Data Release 4 

In this paper we use galaxy samples constructed from 
the New York University Value-Added Galaxy Catalog 
(NYU-VAGC, see Blanton et al. 2005). The version of 
the NYU-VAGC used here is based on the SDSS Data 
Release 4 (Adelman-McCarthy et al. 2006 ) which con- 
sists of three regions: two in the Northern Galactic Cap 
(NGC) and one in the Southern Galactic Cap (SGC). We 
refer the two NGC regions separately as NGCE (which 
is along the celestial equator) and NGCO (which is off 
the equator). The SGC contains three stripes: one along 
the celestial equator, one to the north and one to the 
south. In this paper, we measure the correlation func- 



tions separately for NGCE, NGCO and SGC, and ob- 
tain the variance from these samples. While the total 
correlation functions of the whole samples are obtained 
by combining the pair counts in all these three regions. 
Note that these three regions have quite large angular 
separations which ensure us to obtain a good estima- 
tion of the cosmic variance, except for those faintest 
galaxies at very low redshift where the variance may 
be somewhat underestimated. Our main sample con- 
tains 284489 galaxies with extinction-corrected Petrosian 
magnitudes in the range 14.5 < r < 17.6, redshifts in 
the range 0.01 < z < 0.3 and absolute magnitudes in 
the range -23.0 < M r .o.i < -18.0. The r-band abso- 
lute magnitudes are corrected to redshift z = 0.1 using 
the iVcorrection code of Blanton et al. (2003a) and the 
luminosity-evolution model (.E-correction) of Blanton et 
al. (2003b), given by 

M rfiA = M r , z + 1.62(2 - 0.1) . (1) 

In what follows, all absolute magnitudes and luminosi- 
ties are quoted with such K- and E-corrections. In the 
apparent-magnitude limits adopted here, 14.5 < r < 
17.6, the bright end is chosen to avoid incompleteness due 
to the large angular sizes of galaxies, while the faint end 
is chosen to match the magnitude limit of the SDSS main 
galaxy sample. In our analysis, we will also use another 
apparent-magnitude criterion, 13.0 < r < 17.6 to con- 
struct a supplemental sample. This sample has in total 
287728 galaxies, about 3239 more than the main sample 
described above. Most of these additional galaxies are 
intrinsically bright galaxies with low redshifts, and the 
use of these galaxies can significantly improve the signal 
of the cross-correlation function between the brightest 
and faintest samples by increasing their overlapping vol- 
ume between them. As we will see below, since we can 
normalize the cross-correlation function with the random 
sample constructed according to the faint sample that is 
complete, and since the incompleteness introduced by in- 
cluding the brightest galaxies is not expected to be corre- 
lated with the real structure, the incompleteness should 
not introduce significant bias in our results. 

From our main sample, we first construct 6 volume- 
limited samples (Samples VI - V6) in six bins of 
absolute-magnitude (M r0 i): VI (-23.0,-21.5], V2 (- 
21.5,-21.0], V3 (-21.0,-20.5], V4 (-20.5,-20.0], V5 (-20.0,- 
19.0], V6 (-19.0,-18.0], with detailed selection criteria 
given in Table [T] Obviously, the redshift overlap be- 
tween the brightest and faintest samples is small (and 
sometimes none). In order to enlarge the redshift over- 
lap, which is necessary in order to measure the cross 
correlations reliably, we also construct 6 volume-limited 
samples, Va - Vf, from our supplemental sample. In the 
upper panels of Fig. [T] we show the magnitude-redshift 
distributions of 10% of all these volume-limited sam- 
ples: V1-V6 (upper left panel) and Va - Vf (upper right 
panel) . 

In addition to the volume-limited samples described 
above, we also construct 6 flux-limited samples, Fl - F6. 
As for samples VI - V6, these are defined by the ap- 
parent magnitude limits 14.5 < r < 17.6 and by a bin 
in absolute magnitude. However, no limits are imposed 
on the redshifts, as illustrated in the lower-left panel of 
Fig. [2 The details of these samples are given in Ta- 
ble [2l where the last column indicates the minimum and 
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Fig. 1. — The magnitude-redshift distributions of 10% galaxies in the volume-limited samples VI— V6 (upper left panel), Va - Vf (upper 
right panel) and flux-limited samples Fl - F6 (lower left panel). The lower right panel shows the color-magnitude distribution of 10% of 
the main sample galaxies. The solid line g — r = 0.7 — 0.032(M r .o.i + 16.5) separates the galaxies into red and blue populations. See Tables 
[T]and [2] for more infomation. 



TABLE 1 

Volume-limited samples 



Number of Galaxies 



Sample 




NGCE 


NGCO 


SGC 


z 


Vl/Va 
V2/Vb 
V3/Vc 
V4/Vd 
V5/Ve 
V6/Vf 


(-23.0, -21.5] 
(-21.5, -21.0] 
(-21.0, -20.5] 
(-20.5, -20.0] 
(-20.0, -19.0] 
(-19.0, -18.0] 


4670/5365 
9288/9485 
11833/12110 
10842/11158 
6493/6974 
2620/2690 


8935/10030 
16978/17442 
19987/20542 
17705/18348 
15462/16073 
5321/5485 


2704/3050 
4777/4919 
5676/5837 
5146/5317 
4317/4571 
1534/1590 


(0.098, 0.193)/(0.051, 0.193) 
(0.051, 0.157)/(0.026, 0.157) 
(0.041, 0.127)/(0.021, 0.127) 
(0.033, 0.103)/(0.017, 0.103) 
(0.026, 0.066)/(0.013, 0.066) 
(0.017, 0.043)/(0.010, 0.043) 



Note. — Column 1 indicates the sample ID. Column 2 gives the absolute-magnitude range covered by each sample. Columns 3, 4 and 
5 list the galaxy number in each sample in the three regions (NGCE, NGCO, SGC) of the SDSS. Column 6 lists the redshift range of each 
sample. Each sample is a volume-limited in the redshift range it covers. Samples VI - V6 are for galaxies with r-band apparent magnitudes 
in the range from 14.5 to 17.6. Samples Va - Vf are constructed from all galaxies with r-band apparent magnitudes in the range from 13.0 
to 17.6. The reason for using these two sets of volume-limited samples is given in the main text. 
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TABLE 2 
Flux-limited samples 









Number of Galaxies 




Sample 


M 


NGCE (Red/Blue) 


NGUO (Red/Blue) 


SGC (Red/Blue) z 



Fl (Fl rcd 
F2 (F2 rcd 
F3 (F3 rcd 
F4 (F4 rcd 
F5 (F5 rcd 
F6 (F6 rcd 
p* 



/Flbluc 

/F2 bluc 

/F3biuc 
/F4 bluc 

/F5biu C 
/F6biu C 




8817 (7203/1614) 

14189 (9173/5016) 
19025 (10968/8057) 

17214 (8979/8235) 
20252 (8949/11303) 

5597 (1409/4188) 
18444 



17012 (13745/3267) 
26598 (17150/9448) 
33459 (19174/14285) 
29037 (15186/13851) 
36386 (15810/20576) 
12687 (3570/9117) 
31714 



5108 (4147/961) (0.051, 0.300) 

7556 (4938/2618) (0.041, 0.193) 

9498 (5434/4064) (0.033, 0.157) 

8329 (4271/4058) (0.026, 0.127) 

10162 (4369/5793) (0.017, 0.103) 

3563 (1041/2522) (0.010, 0.066) 
8975 



Note. — Similar to Table [T] but for flux-limited samples. Here again galaxies with r-band apparent magnitudes in the range from 14.5 
to 17.6 are selected. Numbers are listed separately for all, red and blue galaxies. F* is a sample of galaxies with luminosities within a 0.5 
magnitude bin around L* . 
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Fig. 2. — The projected 2PCFs for the flux-limited samples. Here we compare the results obtained by using two different random samples, 
Rl and R2 (see the main text for how Rl and R2 are constructed). The error-bars in this plot represent the l-a scatter among the three 
SDSS regions. For clarity, we only plot error-bars for the auto-correlation functions but note that the error-bars for the cross-correlation 
functions have similar sizes. The results for the cross-correlation functions are shifted downward by 0.5 and 1 dex for the cases Fi — F(i+1) 
and Fi-F(i+2) (i=l,2,3,4,5), respectively. Note that in most of the cases, these two types of random samples yield very similar results. 
However, discrepancy exists in a few cases, especially for the auto-correlation function of F5. As discussed in the main text, this discrepancy 
is likely caused by the supercluter at redshift z ~ 0.08. 
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maximum redshifts of the galaxies in each sample (which 
simply arise because of the apparent magnitude limits). 
The advantage of using flux-limited samples is that we 
have more galaxies in each sample, but the drawback is 
that the samples are not homogeneous in the radial di- 
rection, and so the two samples to be cross-correlated do 
not occupy exactly the same volume. As we will show 
below, the results based on the flux-limited samples are 
very similar to those based on the volume-limited sam- 
ples. Because of this, we will focus on the flux samples 
to take advantage of the larger sample sizes. For each of 
the flux-limited sample, we further divide galaxies into 
two color subsamples according to their g — r colors. The 
dividing line is given by 

g - r = 0.7 - 0.032(M r , .i + 16.5) , (2) 

where g and r are the g-band and r-band magnitudes 
that are if-corrected (Blanton et al. 2003a) and E- 
corrected (Blanton et al. 2003b) to redshift z = 0.1. 
Galaxies with g — r color above (below) this line are 
called red (blue) galaxies, and the corresponding sam- 
ples are denoted by Fl rec i - F6 ro d (red subsamples) and 
Flbiue ~~ F6biuo (blue subsamples). The selection criteria 
for these subsample are listed in Table [U In the lower 
right panel of Fig. Q] we show the color-magnitude rela- 
tion for a random set of 10% of all galaxies in our main 
sample, with the solid line indicating the color separa- 
tion given by Eq. Note that this color separation is 
very similar to what has been used in previous studies 
(e.g. Blanton et al. 2003c; Baldry et al. 2004; Bell et al. 
2004; Hogg et al. 2004; Weinmann et al. 2006). 

As comparison, we also construct a reference sample, 
F*, which includes all galaxies within a 0.5-magnitude 
bin around L*. The selection criteria for this sample are 
listed in the last line of Table [2] 

2.2. The random samples 

In order to measure the two-point correlation func- 
tions, one needs to construct random samples to normal- 
ize the galaxy-pair counts. Since we will study the clus- 
tering properties separately for all, red and blue galaxies, 
we need to generate random samples separately for each 
of them. In general, if we know the luminosity function 
of galaxies, we can first construct a random sample with 
the luminosity function and then apply the observational 
magnitude limits. Unfortunately, we do not have the lu- 
minosity functions for the red and blue galaxies defined 
according to the color criterion described above. Because 
of this, here we adopt a method proposed in Li et al. 
(2006), where a random sample is generated by assign- 
ing each galaxy in the real sample a random position in 
the sky while keeping all other properties (i.e., redshift, 
color, magnitude, etc.) unchanged. In what follows we 
refer to the random samples thus constructed as random 
samples of type Rl. Note that these random samples 
have, by construction, exactly the same redshift distri- 
bution, n(z), as the original sample. Consequently, it is 
not exactly 'random' in that it can still reveal structures 
in n(z). Therefore, this method is only expected to work 
well for wide-angle surveys in which the sky coverage is 
much larger than the large-scale structure in question. In 
addition, it requires that the variation in survey depth 
be small across the sky coverage (Li et al. 2006). 



Although Li et al. (2006) found that two-point corre- 
lation functions based on the random samples described 
above match well those obtained using random samples 
based on the luminosity function, we need to test whether 
or not the use of different random samples can have sig- 
nificant impact on our results for the smaller subsamples 
used here. To do this, we generate another set of random 
samples for all galaxies according to the luminosity func- 
tion obtained by Blanton et al. (2003b). Here a large 
set of points are randomly distributed within the survey 
sky coverage and each point is assigned a redshift and 
a magnitude according to the luminosity function. This 
set of random samples will be referred to as type R2. 

For both sets of random samples, the mangle soft- 
ware provided by A. Hamilton (Hamilton & Tegmark 
2004) is used to locate the polygon ID for each random 
galaxy and to determine the completeness and magni- 
tude limit (which changes slightly across the sky cover- 
age) for this polygon using the mask provided by Blanton 
et al. (2005). The completeness and magnitude limit 
so obtained are then included in the random samples. 
In calculating the two-point correlation functions for all 
galaxies we use random samples that are 10 time bigger 
than the observational data. While for red or blue galax- 
ies we use random samples that are 15 times bigger than 
the observational data. 

3. METHOD TO ESTIMATE THE TWO-POINT 
CORRELATION FUNCTION 

We estimate both the 2-point auto- and cross- 
correlation functions (hereafter 2PCFs) using the defi- 
nition proposed in Davis & Peebles (1983), 

where i, j denote the pair of samples for which the cross 
correlation function is estimated; i = j corresponds to 
auto-correlation function. As a convention, sample i al- 
ways represents the brighter sample if i ^ j. In the above 
definition, DiDj is the count of galaxy pairs between 
samples i and j, DiRj is the count of galaxy-random 
pairs between galaxy sample i and random sample j, and 
r p and rv are the separations of galaxy pairs perpendicu- 
lar and parallel to the line of sight, respectively. Galaxy- 
galaxy and galaxy-random pairs are counted in logarith- 
mic bins in r p , with bin width A\og w (r p ) = 0.2, and in 
linear bins of r w , with bin width Arv = 1 /i Mpc. Since 
we have galaxy samples in three different regions (NGCE, 
NGCO and SGC), we combine the number counts of all 
regions when estimating the correlation function. The 
scatter among the three regions are used to indicate the 
uncertainties (error-bars) in the measurements. Note 
that the errors of the 2PCFs for the faintest sample 
may be somewhat underestimated compared to the true 
cosmic variance because the volumes and separations of 
these three regions are too small to represent the true 
cosmic variance. 

As mentioned in Section 12.11 the redshift overlap be- 
tween some of the volume-limited samples is very small, 
especially between the bright and faint samples. In fact, 
samples VI and V5 as well as VI and V6 have no overlap- 
volume at all (cf. Fig[TJ) . In those cases we have to use 
samples Va - Vf (which may be incomplete) to get better 
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TABLE 3 

The correlation lengths and slopes for volume-limited samples on large scales (0.98 h _1 Mpc < r p < 9.6 /i _1 Mpc) . 



Sample 


Va 


Vb 


Vc 


Vd 


Ve 


Vf 




a 


a 


a 


a 


a 


a 




ro 


ro 


ro 


ro 


ro 


ro 


VI 
V 1 


i qa 4- n f\t\ 
l.o4 it U.UD 














7 £Q 4- fl 7^ 












V 1 


1 SQ 4- n Oft 

l.oo it U.UD 


1 77 4- n n/i 
1. it it U.U4 












ft on 4- n e:n 
D.yU ± u.oy 


D.ll it U.oo 










V 6 


J-.CSo It U.U 1 


1 70 4- n nfi 
1. ( y it U.UD 


l.ol It U.U4 










ft fin 4- n ft/i 
D.DU ± U.U4 


o.oO it U.o4 


K fiO 4- n 1 ft 
O.Dz ± U. ID 








V 1 


l.oz ± U.ll 


1 q 1 _l_ n no 
i.ol it U.Uo 


i ci 4_ n nft 
l.ol ± U.UD 


1 on 4_ n nc; 
l.oU ± U.Uo 








6.41 it 0.51 


5.67 ± 0.37 


5.56 ± 0.25 


5.38 ± 0.16 






V5 


1.71 ± 0.23 


1.79 ±0.16 


1.84 ±0.13 


1.79 ±0.11 


1.83 ±0.08 






6.54 ± 0.94 


5.23 ± 0.64 


5.16 ± 0.40 


4.96 ± 0.32 


4.90 ±0.18 




V6 










1.86 ±0.29 


1.92 ±0.30 












4.70 ±0.43 


4.17 ±0.23 


a = a = 1.82 


ro 


ro 


ro 


^0 


ro 


ro 


VI 


7.57 ± 0.45 












V2 


6.89 ± 0.34 


6.13 ±0.18 










V3 


6.56 ±0.38 


5.87 ±0.20 


5.62 ± 0.11 








VI 


6.40 ± 0.39 


5.67 ±0.23 


5.56 ±0.14 


5.39 ±0.11 






V5 


6.56 ±0.58 


5.26 ±0.37 


5.14 ±0.29 


4.97 ±0.21 


4.89 ±0.09 




V6 










4.65 ±0.38 


4.08 ± 0.20 



Note. — The upper part lists the best fit values of ro and a for different volume-limited samples, using data points on large scales, 
0.98 h -1 Mpc < r p < 9.6 h~ 1 Mpc. The table is not filled up because the overlaps in volume between the faintest and bright samples are 
too small to allow reliable estimates of the cross correlation functions. The lower part of the table shows the correlation length obtained 
by setting the slope of the correlation functions, a, to be a = 1.82, the average of the values listed in the upper part. 



TABLE 4 

The correlation lengths and slopes for volume-limited samples on small scales (0.16 h~ 1 Mpc < r p < 0.98 h _1 Mpc) 



Sample 


Va 


Vb 


Vc 


Vd 


Ve 


Vf 




a 


a 


a 


a 


a 


a 




ro 


ro 


ro 


ro 


ro 


ro 


VI 


1.89 ± 0.16 














8.07 ±0.81 












V2 


1.84 ±0.08 


1.86 ±0.10 












7.82 ± 0.59 


5.77 ±0.49 










V3 


1.83 ±0.08 


1.83 ±0.06 


1.77 ±0.03 










7.89 ±0.66 


5.67 ±0.27 


5.71 ±0.36 








V4 


1.82 ± 0.08 


1.85 ±0.08 


1.71 ±0.03 


1.71 ±0.02 








7.93 ± 0.47 


5.61 ±0.60 


5.88 ±0.56 


5.70 ±0.38 






V5 


1.88 ±0.21 


1.96 ±0.13 


1.81 ±0.09 


1.77 ±0.07 


1.72 ±0.02 






7.62 ± 0.68 


4.90 ± 0.29 


5.27 ±0.38 


5.07 ±0.28 


5.27 ±0.19 




V6 










1.80 ±0.14 


1.75 ±0.13 












4.73 ±0.36 


4.43 ± 0.29 


a = a = 1.81 


r 


r 


ro 


ro 


ro 


ro 


VI 


8.16 ± 0.54 












V2 


7.21 ± 0.40 


6.00 ±0.13 










V3 


7.07 ± 0.36 


5.74 ±0.12 


5.53 ± 0.09 








VI 


6.99 ± 0.46 


5.76 ±0.13 


5.44 ± 0.08 


5.27 ±0.06 






V5 


6.70 ±0.26 


5.45 ±0.23 


5.25 ± 0.12 


4.93 ± 0.09 


4.90 ±0.11 




V6 










4.69 ±0.14 


4.23 ± 0.09 



Note. — The same as Table(3] except that the results are obtained by fitting the correlation functions on small scales, 0.16 h x Mpc < 
r p < 0.98 h~ 1 Mpc. 



statistics. According to Eq. ([3]), one may circumvent the 
problems to some extent using combinations of complete 
and incomplete samples. For example, one can choose 
the i sample from Va - Vf, and the j sample and the 
corresponding random sample from VI - V6. Since VI 
- V6 are complete, this reduces the effects due to the 
incompleteness in samples Va - Vf. 
In our following discussion, we will focus on the pro- 



jected 2PCF, which is defined as 



Jo k 

(4) 

Here the summation is made over k — 1 to 40 corre- 
sponding to rv = 0.5 /i _1 Mpc to rv = 39.5 /i _1 Mpc. If 
we assume that the real-space 2PCF £(r) can be fitted 
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Different Wy 

Fig. 3. — The correlation length, ro, and logarithmic slope, a, 
obtained from the projected auto- and cross-correlation functions 
on scales 0.98 h~ 1 Mpc < r v < 9.6 h~ x Mpc for various volume- 
limited and flux-limited samples. The error-bars are 1-cr scatter 
among the three SDSS regions. In the horizontal axis we indicate 
the IDs (i — j) of the samples used in the estimate of the corre- 
lation function, Wij. The circles with error-bars are the results 
for volume-limited samples. The squares with error-bars are the 
results for flux-limited samples. Note that since the overlap region 
is too small, we do not have measurement of 1-6, 2-6, 3-6 and 4- 
6 data points for volume-limited samples and 1-6 for flux-limited 
samples. Overall, the results for volume-limited and flux-limited 
samples agree well. Large discrepancy exists for Wg g, likely due 
to the existence of the supercluster at redshift z ~ 0.08. 



by a power law, 



iir) = (r„/r) c 



(5) 



where r$ and a are, respectively, the correlation length 
and slope of the correlation function, then the projected 
2PCF is related to the real-space 2PCF as 



W(r p ) 



'O'p 



, a r (l/2)r[(q-l)/2] 
r(a/2) 



(6) 



where T(x) is the Gamma function. As we will see below, 
the observed 2PCFs can all be fitted reasonably well by 
power laws over limited ranges of r p . We can then use ro 
to represent the correlation amplitude, and a the shape 
of the correlation function. 

4. DEPENDENCE ON GALAXY LUMINOSITY 

We start our investigation by comparing the pro- 
jected 2PCFs obtained using the volume-limited and 
flux-limited samples. When calculating the projected 
2PCFs with the volume-limited samples, we use com- 
binations of the set V1-V6 with the set Va-Vf for the 
reasons given in Section [3J Since the results are very sim- 
ilar for the volume-limited and the flux-limited samples, 
we show only the results for the flux-limited samples in 
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Fig. 4. — The same as Fig. \3\ except that the results are ob- 
tained from the fitting to the correlation functions on small scales 
(0.16 /i~ 1 Mpc < r p < 0.98 /i~ 1 Mpc). Note again that since the 
overlap region is too small, we do not have measurement of 1-6, 
2-6, 3-6 and 4-6 data points for volume-limited samples and 1-6 for 
flux-limited samples. 



Fig. [21 where we compare the projected 2PCFs based on 
random samples of type Rl with those based on random 
samples of type R2. As found by Li et al. (2006), the two 
types of random samples give almost identical results in 
most cases. The largest discrepancy is seen in the auto 
2PCF of sample F5, This sample is significantly affected 
by the supercluter at redshift z ~ 0.08, and the points in 
sample Rl is not completely random due to the way it 
is constructed. Since overall the random samples of type 
Rl work quite well and since their construction does not 
require the luminosity function of galaxies, the results 
presented below will be based on this type of random 
samples. A comparison of the various projected 2PCFs 
in Fig.[2]shows that brighter galaxies are on average more 
strongly clustered, an effect that we will describe in more 
detail in the following. 

To quantify the observed 2PCFs, we fit to the projected 
2PCFs according to Eq. © with the two free parame- 
ters ro and a. A larger value of ro implies a stronger 
clustering strength, while a larger a implies a steeper 
correlation function. In order to see how ro and a might 
change with scales, we fit the projected 2PCFs separately 
over two scale ranges, 0.98 ft, _1 Mpc < r p < 9.6 /i _1 Mpc 
and 0.16 /i^Mpc < r p < 0.98 /i _1 Mpc. The separation 
between large and small scales is made at r p ~ l/i~ 1 Mpc, 
because the correlation function on smaller scales is ex- 
pected to be dominated by the one-halo term, while that 
on larger scales by the two-halo term, as shown in Yang 
et al. (2005). We use two methods to obtain ro and a 
from the data. In the first method, we obtain ro and a 
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Fig. 5. — The relative bias as a function of r p obtained from the ratios of the projected correlation functions of the flux-limited samples. 
In the lower-right panel, we show the mean relative bias parameters obtained by averaging the results shown in each of the five panels. In 
each panel, i — j means that the projected two-point correlation function is for samples Fi and Fj. 



TABLE 5 

The correlation lengths and slopes for flux-limited samples on large SCALES (0.98 /l -1 Mpc < T p < 9.6 h _1 Mpc) 



Sample 


Fl 


F2 


F3 


F4 


F5 


F6 




a 


Q 


a 


a 


a 


a 




ro 


ro 


ro 


ro 


ro 


ro 


Fl 


1.81 ± 0.03 














7.83 ±0.61 












F2 


1.80 ±0.03 


1.76 ±0.03 












6.88 ± 0.42 


6.06 ±0.22 










F3 


1.84 ±0.04 


1.77 ±0.04 


1.73 ± 0.02 










6.54 ± 0.33 


5.82 ±0.22 


5.91 ± 0.14 








F4 


1.80 ±0.07 


1.77 ±0.05 


1.79 ±0.03 


1.78 ±0.03 








6.36 ± 0.63 


5.73 ±0.27 


5.47 ±0.11 


5.41 ±0.13 






F5 


1.81 ± 0.13 


1.75 ±0.07 


1.77 ±0.05 


1.78 ± 0.03 


1.64 ±0.02 






6.09 ±0.88 


5.46 ±0.31 


5.29 ±0.17 


5.17 ±0.09 


5.55 ± 0.15 




F6 




1.76 ±0.35 


1.86 ±0.20 


1.83 ±0.16 


1.85 ±0.09 


1.88 ±0.12 






5.59 ± 1.12 


5.14 ±0.48 


4.96 ±0.34 


4.73 ± 0.14 


4.49 ±0.13 


a = a = 1.79 


ro 


ro 


ro 




ro 


ro 


Fl 


7.80 ± 0.31 












F2 


6.86 ±0.21 


6.09 ±0.13 










F3 


6.49 ±0.19 


5.83 ±0.12 


5.94 ±0.14 








Fl 


6.34 ±0.35 


5.75 ±0.18 


5.48 ±0.15 


5.42 ±0.14 






F5 


6.07 ±0.56 


5.48 ± 0.29 


5.31 ± 0.10 


5.17 ±0.09 


5.65 ± 0.11 




F6 




5.62 ±0.46 


5.08 ± 0.37 


4.93 ±0.21 


4.68 ± 0.11 


4.38 ±0.12 
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TABLE 6 

The correlation lengths and slopes for flux-limited samples on small scales (0.16 h~ 1 Mpc < r p < 0.98 /i _1 Mpc) 



Sample 


r 1 


T?0 

r A 


r o 


r 4 


r O 


r O 




ft 


Ol 


ft 


ft 


ft 


Ol 




Tn 
' u 


Tn 
' U 


Tn 
' u 


' U 


Tn 
' U 


Tn 
' U 


Fl 


i on _|_ n no 














8.31 ± 0.75 












F2 


1.82 ± 0.04 


1.87 ±0.05 












7.88 ±0.41 


5.72 ± 0.39 










F3 


1.83 ±0.03 


1.82 ±0.04 


1.72 ± 0.04 










7.85 ± 0.34 


5.61 ±0.27 


5.97 ±0.36 








F4 


1.81 ±0.04 


1.83 ± 0.06 


1.72 ± 0.03 


1.72 ±0.02 








7.84 ± 0.41 


5.55 ±0.26 


5.66 ±0.18 


5.62 ±0.17 






F5 


1.81 ±0.04 


1.88 ±0.06 


1.76 ±0.03 


1.72 ± 0.02 


1.69 ± 0.03 






7.66 ±0.51 


5.19 ±0.13 


5.42 ± 0.14 


5.39 ±0.16 


5.68 ±0.19 




F6 




1.96 ± 0.32 


1.83 ± 0.12 


1.80 ±0.08 


1.73 ± 0.03 


1.75 ±0.04 






5.07 ± 1.21 


5.43 ± 0.45 


5.15 ±0.44 


5.11 ± 0.21 


4.82 ±0.20 


a = a = 1.80 


7-0 




ro 




ro 


ro 


Fl 


8.21 ± 0.42 












F2 


7.17 ±0.23 


6.01 ±0.19 










F3 


6.93 ±0.11 


5.71 ±0.19 


5.61 ± 0.15 








F4 


6.96 ± 0.13 


5.69 ±0.18 


5.33 ± 0.06 


5.29 ±0.15 






F5 


6.77 ±0.16 


5.50 ±0.15 


5.25 ± 0.17 


5.10 ±0.18 


5.24 ±0.14 




F6 




5.69 ±0.59 


5.54 ±0.09 


5.13 ±0.13 


4.85 ± 0.17 


4.64 ±0.15 



TABLE 7 

The correlation lengths and slopes for flux-limited samples of red galaxies on large scales 

(0.98 /i^Mpc <r p < 9.6 h^Mpc) 



Sample 


Fired 


F2 rcd 


F3red 


F4 rcd 


F5 ro d 


F6 ro d 




a 


a 


a 


a 


a 


a 




ro 


ro 


ro 


ro 


ro 


ro 


Fired 


1.83 ± 0.03 














8.50 ± 0.59 












F2 rcd 


1.84 ±0.05 


1.80 ±0.04 












7.59 ± 0.48 


6.82 ±0.26 










F3 rcd 


1.89 ±0.05 


1.82 ±0.04 


1.87 ±0.03 










7.33 ± 0.41 


6.70 ±0.24 


6.53 ± 0.20 








F4 rcd 


1.85 ± 0.08 


1.84 ±0.05 


1.87 ±0.04 


1.88 ±0.03 






7.34 ± 0.68 


6.78 ± 0.30 


6.59 ±0.18 


6.71 ±0.18 






F5 rcd 


1.85 ± 0.15 


1.86 ±0.07 


1.90 ±0.05 


1.90 ± 0.03 


1.95 ± 0.02 




7.23 ± 1.03 


6.71 ±0.36 


6.63 ± 0.24 


6.66 ±0.12 


6.55 ± 0.11 




F6 rcd 




1.91 ±0.24 


2.01 ± 0.19 


1.95 ±0.13 


2.08 ±0.07 


2.15 ±0.07 






7.19 ±1.11 


7.16 ± 0.78 


6.93 ±0.56 


6.89 ± 0.18 


7.18 ±0.18 


a = a = 1.90 


ro 


r 


r () 


ro 


ro 


ro 


Fired 


8.50 ±0.44 












F2 rcd 


7.60 ± 0.36 


6.82 ±0.25 










F3 rcd 


7.33 ± 0.35 


6.71 ±0.25 


6.52 ± 0.13 








F4 rcd 


7.35 ± 0.37 


6.77 ±0.34 


6.59 ±0.21 


6.70 ±0.14 






F5 rc d 


7.24 ±0.32 


6.70 ± 0.33 


6.62 ± 0.23 


6.66 ±0.12 


6.56 ± 0.09 




F6 rcd 




7.19 ±0.78 


7.17 ±0.19 


6.92 ±0.22 


6.97 ±0.29 


7.25 ±0.31 
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TABLE 8 

The correlation lengths and slopes for flux-limited samples of red galaxies on small scales 

(0.16 ft-^pc < r p < 0.98 h~ 1 Mpc) 



Sample 


Fl rG d 


F2 ro d 


F3 rc d 


F4 rod 


F5 rc d 


F6 rod 




a 


a 


a 


a 


a 


a 






r 


r 


r 


r 






i qi x n in 














q -1- n on 
o.oy ± u.yu 












r 2 rc d 


i flft 4- n n7 
l.oD ± u.u/ 


l.yo it U.UD 












o.4u i u.oy 


ft Ofl 4- n OQ 

D.zo it U.zo 










T?q 

1 °Tcd 


i e7 _i_ n nft 
i.oi zt U.UO 


l.yz it U.U4 


i qa 4- n no 
J..CS4 It U.Uz 










o.oo ± U.DO 


ft Qf; 4_ n q /i 
O.oD it U.o4 


ft ftO 4_ n OQ 
U.OZ it U.zo 








F4,„i 


1.90 ± 0.08 


1.94 ± 0.06 


1.82 ± 0.03 


1.84 ± 0.02 








8.27 ±0.89 


6.40 ± 0.38 


6.69 ±0.44 


6.94 ±0.27 






F5 rcd 


1.91 ±0.27 


1.86 ±0.22 


1.85 ±0.18 


1.82 ±0.12 


1.85 ±0.07 






8.28 ± 1.12 


8.40 ± 0.34 


7.20 ±0.36 


7.50 ± 0.24 


7.59 ±0.16 




F6 rc d 




1.91 ±0.20 


1.90 ±0.27 


1.86 ±0.29 


1.88 ±0.16 


2.04 ±0.19 






8.17 ± 1.64 


8.66 ±0.66 


8.39 ±0.41 


8.85 ±0.29 


8.58 ±0.26 


a = a = 1.89 




ro 


ro 


^0 


ro 


ro 


llrcd 


8.52 ± 0.49 












F2 rcd 


7.62 ± 0.43 


6.61 ±0.18 










F3 rc d 


7.75 ± 0.33 


6.54 ±0.14 


6.42 ±0.31 








F4 rcd 


8.04 ± 0.23 


6.72 ±0.16 


6.59 ±0.33 


6.76 ±0.26 






F5 rcd 


8.13 ± 0.49 


7.03 ±0.18 


6.97 ±0.16 


7.00 ± 0.24 


7.33 ±0.15 




F6 rcd 




7.63 ±0.83 


8.78 ±0.36 


8.04 ±0.34 


8.64 ±0.22 


10.12 ±0.25 
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Fig. 6. — The relative bias parameters obtained from the cor- 
relation functions over 0.98 h~ 1 Mpc to 9.8 h~ 1 Mpc using seven 
different approaches (see the main text for details). The error- 
bars are based on the scatter of the data points in the range 
0.98 h _1 Mpc <r p < 9.8 /i~ 1 Mpc. Panels from top to bottom cor- 
respond to the bias parameters for Fi (i = 1, 2, 3, 5, 6), normalized 
by that of F4. The horizontal axis labels the different approaches 
that are used to determine the bias parameter for each subsample. 
One data point is missing in the first and fifth panels, because the 
redshift overlap between Fl and F6 is too small to measure the 
cross-correlation function reliably. 



directly from the fit to the projected 2PCFs. In the sec- 
ond one, we first obtain an average value of a, a, from all 
the correlation functions and then fit the projected 2PCF 
to get ro keeping a — a. The fitting parameters so ob- 
tained are listed in Tables [21 HI El and El separately for 
volume-limited and flux-limited samples, and for large 
and small scales. 

Fig. [3] compares the values of ro and a on large scales 
(0.98 /i _1 Mpc < r p < 9.6 ft^Mpc) obtained from the 
volume-limited samples with those obtained from the 
corresponding flux-limited samples. Overall, the results 
for these two kinds of samples are similar. The only ex- 
ception is PFs^, where the values of ro and a for the 
volume-limited and flux-limited samples are significantly 
different, largely because the presence of the supercluster 
at z ~ 0.08 that makes the correlation properties quite 
sensitive to how the sample is exactly formed. There is 
a clear trend in the correlation length ro that brighter 
galaxies have larger ro. The cross-correlations between 
bright and faint galaxies are also stronger than the auto 
correlations of faint galaxies. On the other hand the 
slope of the correlation function does not show any strong 
trend with luminosity, and the mean value of and scatter 
in a are 1.82 ± 0.04 for the volume-limited samples and 
1.79 ± 0.05 for the flux-limited samples. The correlation 
lengths obtained by fixing a at these mean values are 
also listed in Tables El [4] El and El for comparison. 

In order to investigate how the clustering strengths 
change on small scales, we fit the projected 2PCFs 
with power laws in the range 0.16 /i _1 Mpc < r p < 
0.98 h~ 1 Mpc, and the results are shown in Fig. [4] and 
listed in Tables [4] and El A comparison to the results 
shown in Fig. E] shows that the clustering strengths for 
the samples that involve the brightest galaxies are en- 
hanced on small scales; in particular the cross correla- 
tion between the brightest and faintest galaxies is much 
stronger on small scales. For faint galaxies, the clustering 
strengths on small scales are similar to those based on the 
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Fig. 7. — Left panel: The average relative bias parameter, obtained using the correlation data in 0.98 h _1 Mpc < r p < 9.8 h~ 1 Mpc as 
a function of galaxy luminosity. The results are obtained using the flux-limited samples, and both the bias parameter and the luminosity 
are normalized by the values for F4. For a given luminosity, the mean and error-bar are based on the different measurements shown in 
Fig|6] The solid curve shows the relation, bi/64 = 0. 862 + 0. 120L{ / L4, obtained by fitting the data. Right panel: In this panel, we compare 
our results with those obtained by Zehavi et al. (2005). Here the bias parameters and the luminosities are normalized by the values for 
L* galaxies. In Zehavi et al., the relative bias parameters are measured using the auio-correlation functions at r p = 2.7 h _1 Mpc for six 
volume-limited samples with luminosities in (-23.0, -22.0), (-22.0, -21.0), (-21.0, -20.0), (-20.0, -19.0), (-19.0, -18.0), (-18.0, -17.0), 
and their results are shown by triangles. Using the same method and same samples, our results (however at r p = 2.8 h —1 Mpc) are shown 
as open circles and match very well their results. The open squares show the same results as shown in the left panel, except that the 
normalizations are different. The curve is a fit to the results represented by the open squares. 



power laws obtained on large scales. This suggests that 
there is a significant population of faint galaxies that are 
distributed around the brightest galaxies and the num- 
ber drops rapidly with the increase of the distance to 
the brightest galaxies. As we will see below, it is likely 
that this population is dominated by faint red galaxies 
in clusters and groups. 

Having investigated the clustering strengths for galax- 
ies of different luminosities, we now study the bias factor 
in the distribution of galaxies as a function of galaxy lu- 
minosity. Since we do not know the two-point correlation 
function of the cosmic density field, we can only study 
the bias factor in a relative sense. To do this, we choose 
F4 as our fiducial sample, because it has a good overlap 
in volume with both the faint and bright samples. Since 
each sample is used in a number of projected 2PCFs, we 
can estimate the relative bias for each sample in a num- 
ber of ways. With F4 as the reference sample, the bias 
factor for the ith sample can be estimated from one of 
the following definitions: 

£ = ^0 = 1,2,3,4,5,6), (7) 

where i = 1,2, 3, 5, 6. Thus, each relative bias can be esti- 
mated using seven different approaches. In what follows, 
these seven approaches will be referred as approaches I, 
II, • • •, VI, corresponding to j = 1, 2, • • •, 6, respectively, 
and approach VII, corresponding to the definition in the 
second line of equation ([7]). Note that in general the 
relative bias factors, &i/&4, obtained from different ap- 
proaches are not expected to be the same. As we will 
discuss in Section [6j these relative bias factors are ex- 



pected to be the same only under the assumption that 
the bias for galaxies of different luminosities is linear and 
that the stochastic component in the bias relations be- 
tween different galaxies is small. 

The first 5 panels in Fig. \5\ show the relative bias 
factor as a function of r p for the five samples, Fi 
(i = 1,2,3,5,6), and the 7 different lines within each 
panel show the results obtained using the 7 different ap- 
proaches described above. For most cases, the bias func- 
tion is almost independent of r p at r p > 1 /i _1 Mpc, in- 
dicating that nonlinearity in the relative bias is small on 
these scales. The bias factors obtained using the 7 differ- 
ent approaches agree quite well with each other. Their 
averages are shown in the lower-right panel which clearly 
shows that brighter galaxies have higher relative biases. 

To quantify the relative bias as a function of luminos- 
ity, we estimate a mean relative bias factor for each case 
by averaging the bias function over l/i _1 Mpc-10/i~ 1 Mpc. 
The results are shown in Fig. [5] for different luminosity 
samples using the 7 approaches. Note that there is one 
data point missing in the first and fifth panels, because 
the overlap between samples Fl and F6 is too small to 
have a reliable measurement of Wi t Q. The dashed line in 
each panel is the average value of the relative bias ob- 
tained from the seven different approaches. As one can 
see, the relative bias factors obtained with different ap- 
proaches are consistent with each other. We estimate the 
mean luminosity of galaxies in each sample F1-F6 as 
and the corresponding average absolute magnitudes thus 
obtained are -21.87, -21.24, -20.76, -20.27, -19.61, 
and —18.62, respectively. These magnitudes are then 
converted into luminosities to represent the characteris- 
tic luminosities of the samples in question. The left panel 
of Fig. [7] shows the average relative bias factor versus 
Li/Li as open circles, with error-bars representing the 
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Fig. 8. — Here we compare the correlation lengths (ro) and slopes (a) of all (circles), red (squares) and blue (triangles) galaxies using 
flux-limited samples (Ff - F6 for all galaxies; Fl rct i - F6 re< j for red galaxies, and Fl^iue — F6bi ue f° r blue galaxies). Here the results are 
based on correlation functions on large scales, 0.98 h~ x Mpc < r p < 9.6 h~ 1 Mpc. To avoid confusion, data points for rad and blue galaxies 
are shifted slightly horizontally. 



scatter among the seven different approaches. There is 
a clear trend that the relative bias increases with galaxy 
luminosity. The solid curve is a linear fit to the data: 

bi/bi = 0.862 + 0.120L l / J L 4 - 

The luminosity dependence of galaxy bias in the 
SDSS has been investigated by Zehavi et al. (2005) 
using the amplitudes of the auto correlation func- 
tions at a fixed projected radius r p = 2.7 h~ 1 Mpc 
obtained from six volume-limited samples with M r 
in (-23,-22), (-22,-21), (-21,-20), (-20,-19), 
(—19,-18), (—18,-17). The bias factors were normal- 
ized relative to L* galaxies defined to be the ones with 
M r in (—21, —20). To compare our results with theirs, we 
show in the right panel of Fig. [7] the relative bias factors 
obtained by Zehavi et al. as triangles together with our 
results (open circles) that are obtained from the auto- 
correlation amplitudes at a fixed, but a slightly different, 
projected radius r p = 2.8 ft, _1 Mpc. Note that our results 



are in excellent agreement with the results obtained by 
Zehavi et al. (2005). For the relative bias - luminosity 
relation shown in the left panel of Fig. [7| we can convert 
it into a relation with L/L* . The characteristic luminos- 
ity for sample F* is —20.47, very close to L* . Using the 
projected 2PCF for F*, we can write 



6(/i , = i)(/ ^_M, <„ 

where and W* t * are the projected auto 2PCFs of 

samples F4 and F*, respectively. We first calculate 
bi/b* at different scales and then average these biases 
over 1 ft^Mpc-lO /i~ 1 Mpc to obtain the average bias 
bi/b*. The resulting bi/b* - L/L+ relation is shown by 
the open squares in the right panel of Fig. [7] Fitted 
with a linear model, this relation can be described as 
bi/b* = 0.872 + 0.134Lj/L*. As one can see, our results 
based on the cross-correlation match well those of Zehavi 
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Fig. 9. — The same as Fig. [8] except that the results are based 
0.98?i- 1 Mpc. 

et al. (2005). The implications of this agreement will be 
discussed in Section [5] 

5. DEPENDENCE ON GALAXY COLOR 

In order to study how galaxy clustering depends on 
both galaxy luminosity and color, we have subdivided 
each flux-limited sample, Fi, into a red and a blue sub- 
sample according to the criteria described in equation 
(|5J). We remind that these subsamples are referred to 
as Fi re d and Fibi ue - We have carried out the same anal- 
yses as described in Section |4] for these color samples. 
Figs. [5] shows the correlation lengths ro and slopes a ob- 
tained by fitting the correlation functions in the range 
0.98 /i _1 Mpc < r p < 9.6 /i~ 1 Mpc; the corresponding re- 
sults for small scales, 0.16 h^Mpc < r p < 0.98 h^Mpc, 
are shown in Figs. O In both figures, we also include the 
results for the total samples (i.e. without color separa- 
tion) for comparison. Again, the error-bars are obtained 
from the scatter among the three regions in the SDSS ob- 
servations. To make the results accessible to others, we 



on the correlation functions on small scales, 0.16 h J Mpc < r p < 

list all the values of ro and a in Tables [7] (red galaxies; 
large scales), |5] (red galaxies; small scales), M (blue galax- 
ies; large scales), and [TOl (blue galaxies; small scales). In 
all these tables, we also include the best fit values for ro 
by fixing a — a. 

As one can see clearly from the left panels in Fig. [51 on 
large scales red galaxies in all cases have larger ro than 
all galaxies and blue galaxies in the same luminosity bin. 
Note also that the luminosity dependence of ro is seen 
in all, red and blue samples, although the dependence 
is stronger for blue galaxies, as to be quantified in the 
following. Overall, the correlation slope, a (shown in the 
right panel of Fig. [5]), is larger for red galaxies, and the 
dependence is stronger for the cases where at least one 
of the two samples in cross-correlation is faint. 

As shown in Fig. [51 the correlation lengths ro obtained 
from the projected 2PCFs on small scales have properties 
similar to those on large scales, in that red galaxies have a 
larger correlation length. But the trend is quite different 
in detail. Firstly, the difference between the correlation 
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Fig. 10. — The correlation lengths (ro) and slopes (a) of the cross-correlation functions between red and blue galaxies (squares) on 
large scales, 0.98 h~ 1 Mpc < r p < 9.6 h~ x Mpc. Panels from top to bottom correspond to flux limited samples for blue galaxies Fit,i ue 
(i = 1, 2..., 6). The red samples used in the cross-correlations are labeled on the horizontal axis Fi re( j (i = 1, 2..., 6). For comparison, results 
for all galaxies in the corresponding flux-limited samples are plotted as open circles. 



lengths obtained from the cross correlations of the bright- 
est red sample (Fl re( j) with other red samples are signif- 
icantly larger here than that shown in Fig. [5] Since the 
brightest red galaxies are located predominately in clus- 
ters and rich groups, this suggests that many of the faint 
red galaxies are also in clusters and groups. Secondly, 
for F3 rc d through F5 re d, their cross-correlation with the 
faintest red galaxies has an amplitude that is significantly 
larger than their cross-correlation with brighter (except 
the brightest) red galaxies. This indicates that many 
red galaxies of intermediate luminosities may possess ha- 
los of red satellite galaxies that are more than 1 magni- 
tudes fainter. Finally, on small scales, the luminosity- 
dependence ro for blue galaxies appears to be weaker. 

It is also interesting to look at the cross correlation be- 
tween red and blue galaxies. In Figs.flOland llll we com- 
pare the values of ro and a for such cross correlations on 
large and small scales with those for all galaxies. These 



values are also listed in Tables [TT] and [T2j As one can 
see, on large scales (0.98 /i -1 Mpc < r p < 9.6 /i -1 Mpc) 
(Fig. [TO]), the values of r and a for the red/blue cross- 
correlation functions are quite similar to those for all 
galaxies. However, there is a marked difference on small 
scales (O.ie/i^Mpc <r p < 0.98/i _1 Mpc) (Fig.HJ). Here 
the cross correlation amplitudes between the brightest 
red galaxies (Fl re( j) and blue galaxies are significantly 
lower than those for all galaxies. Since many of the 
brightest red galaxies are the brightest central galaxies 
in clusters and groups, this result suggests the satellite 
galaxies of these galaxies are preferentially red, in qual- 
itative agreement with the findings by Wcinmann et al. 
(2006). 

In order to quantify how galaxies of different colors 
and luminosities are biased with respect to each other 
on large scales, we measure the relative bias factors for 
red and blue galaxies based on the ratios of the corre- 
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Fig. 11. — The same as Fig. 1101 but here results are based on correlation functions on small scales, 0.16 /i _1 Mpc < r p < 0.98 h~ 1 Mpc. 



sponding projected 2PCFs, as we did in Section [4] for 
the luminosity samples. Note that, in addition to the 
red-red and blue-blue cross-correlation functions among 
the 6 luminosity bins, one can also use the red-blue cross- 
correlation functions to obtain the relative bias. These 



are 



W, 



Ri.Bj 



for red galaxies, and 



Wn4,Bj 



W B j,Rj 
WB4,Rj 



(9) 



(10) 



for blue galaxies, where i = 1,2,3,5,6 and j = 
1, 2, 3, 4, 5, 6, and Bi and Ri denotes Fibiue and Fi rcc i, re- 
spectively. In what follows, these six approaches will be 
referred as approaches VIII, IX, ■ • ■, XIII, corresponding 
to j = 1, 2, 6, respectively. Together with the previous 
seven approaches, we have now in total 13 approaches 
to measure the relative bias bi/b^ for each red or blue 



subsamples. Figs. [12] and [13] show the large-scale rela- 
tive bias factors obtained with the 13 approaches for red 
and blue galaxies, respectively. Some points are missing 
again because the overlaps of the samples in question are 
too small. The dashed line in each panel is the average 
value of the relative bias. As one can see, most of the ap- 
proaches give extremely consistent results. The biggest 
outlier is the one that corresponds to the cross correla- 
tion F2 rc( j-F6biuc (approach XIII in the second row of 
Fig I12|) . which deviates from the mean bias ratio by just 
over 2a. 

In the left panels of Fig. [14] (upper panel for red 
galaxies and lower panel for blue galaxies), we show the 
average relative bias (6^/64 obtained from the 13 ap- 
proaches) versus luminosity L2/L4, where Lj for each 
sample is estimated in the same way as described in 
Section [4] Fitting to the data with a linear model 
gives 6i/i»4 = 0.953 + 0.057Li/ L4 for red galaxies, and 
bi/bi = 0. 860 + 0. 132Li/L4 for blue galaxies. As one can 
see, the luminosity-dependence is much stronger for blue 
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TABLE 9 

The correlation length and slope for flux-limited samples of blue galaxies on large scales 

(0.98 h~ 1 Mpc <r p < 9.6 ft.- 1 Mpc) 



Sample 




Flblue 
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F5blue 
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r Oblue 
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A flu 4_ n Aft 
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i en 4- n no 










4.97 ±0.45 


4.62 ±0.21 


4.39 ± 0.32 


4.14 ±0.43 






F5blue 




1.78 ±0.18 


1.65 ±0.13 


1.69 ±0.09 


1.67 ±0.12 


1.67 ±0.15 








4.85 ±0.41 


4.29 ± 0.46 


4.13 ±0.27 


3.90 ±0.38 


3.71 ±0.44 




F6blue 






1.64 ± 0.04 


1.78 ±0.04 


1.77 ±0.04 


1.73 ±0.04 


1.74 ±0.09 








4.71 ± 1.49 


4.10 ±0.57 


3.93 ±0.35 


3.61 ± 0.34 


3.46 ± 0.33 


a = a = 


1.74 


ro 


ro 


ro 


ro 


ro 


ro 


Flblue 




6.15 ± 0.58 












F2 bluc 




5.47 ±0.23 


5.08 ± 0.26 










F3biu C 




4.99 ±0.17 


4.75 ±0.21 


4.62 ± 0.24 








F4 bluo 




4.81 ±0.21 


4.66 ±0.31 


4.42 ± 0.24 


4.22 ±0.23 






F5bi uc 




4.77 ±0.36 


4.47 ±0.29 


4.21 ±0.21 


4.03 ±0.18 


3.80 ± 0.13 




F6biuc 






4.90 ±0.69 


4.02 ± 0.23 


3.89 ±0.37 


3.62 ± 0.19 


3.4b ± 0.13 










TABLE 10 












The correlation length and slope for flux-limited samples of blue galaxies on small scales 










(0.16 h~ 


iMpc < r p < 0.98 h~ 
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1.81 ±0.11 


1.64 ±0.06 


1.65 ±0.06 










5.30 ± 0.45 


4.14 ±0.43 


4.21 ±0.23 


4.06 ±0.38 






F5blue 




1.99 ±0.33 


1.81 ±0.14 


1.64 ±0.07 


1.61 ±0.08 


1.62 ± 0.07 








4.11 ±0.88 


3.86 ± 0.49 


4.03 ± 0.41 


3.88 ±0.38 


3.60 ±0.39 




F6blue 






1.62 ±0.26 


1.73 ±0.32 


1.71 ±0.21 


1.66 ±0.30 


1.60 ±0.28 








5.10 ± 1.92 


3.98 ± 0.97 


3.91 ±0.52 


3.47 ±0.42 


3.45 ±0.51 


a = a = 


1.70 


r 


ro 


ro 


ro 


ro 


ro 


Flbluc 




6.13 ± 0.96 












F2biu C 




5.96 ± 0.26 


5.22 ±0.25 










F3bi uc 




5.23 ±0.38 


4.71 ±0.27 


4.42 ± 0.28 








F4 b i uc 




5.33 ±0.35 


4.40 ± 0.22 


4.09 ±0.16 


3.96 ±0.22 






F5bi uo 




4.77 ±0.35 


4.10 ±0.27 


3.92 ± 0.24 


3.73 ±0.23 


3.53 ± 0.15 




F6bi uc 






4.89 ±0.66 


4.05 ± 0.22 


3.94 ±0.17 


3.44 ±0.17 


3.33 ±0.13 



galaxies than for red galaxies. For red galaxies, signifi- 
cant luminosity-dependence is seen only at the brightest 
end, and there is indication that the bias factor becomes 
larger for the faintest sample. 

In order to put the bias factors for different galaxies 
in the same absolute scale, we normalize all bias factors 
relative to that for all the galaxies (blue and red) in the 
F* sample, i.e. we define 



h/b* = h/bJ 



for red galaxies, and 

for blue galaxies. Here W^^a, Wba,ba and W*'* are the 
projected auto 2PCFs for FA re d, FAuue and F*, respec- 
tively. The values of bi/b* versus L/L* are shown in the 
right panels of Fig. [14] (upper panel for red galaxies and 
lower panel for blue galaxies). Normalized in this way, 
the bias factors for red galaxies are all larger than 1, while 
those for blue galaxies smaller than 1, for the luminosity 
range probed here. Fitting the luminosity dependence 
with a linear model, we get bi/b* = 1.352 + 0.096Li/i* 



The cross-correlation between galaxies 
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TABLE 11 

The correlation length and slope between red and blue galaxies on large scales (0.98 h~ 1 Mpc < r p < 9.6 h~ 1 Mpc) 



Sample 




Fired 


F2 re d 


F3 rcc i 


F4 rc d 


F5 re d 


F6 re a 






a 


Q 


a 


a 


a 


a 






ro 


ro 


ro 


ro 


ro 




* Iblue 




i 7Q -1- n r\K 
1. /y ± U.Uo 


1 ci _i_ n t\A 
l.ol It U.U4 


i en -U n D/i 

l.oU it U.U4 


1 SO 4- n (iA 

i.oy it U.U4 


i oi 4- n 1 n 
l.yi it U.1U 








u.ou zn u.y4 




X. QQ _L_ fl AA 


p. i p. 4_ n ^ 

U. ID It U.OU 


ft. 1 fi 4- 1 A A 
O . lO ZH 1 . 44 




r ^bluo 




1. (4 it U.Uo 


1 71 4.(1 no 

l.i 1 ± u.uz 


i 7ft -U n no 
1. 1 rt U.Uz 


1 7Q 4- n nQ 
l./o ± U.Uo 


i 71 4- n nft 
1. / 1 It U.UD 


1 ft/i _i_ n i q 
1.D4 it U.lo 






ft on r. 4-i 
O.ZD ± U.41 


k fti i n 1B 

0.01 it U.lO 


0.40 ± u.iy 


a q 4_ n 07 
0.4o it U.z / 


K Q/1 4_ n A K 
0.O4 ± U.4D 


O.U1 ± 1. (U 


T?Q 

r Oblue 




1 70 -1- fl HQ 

1. 1 y it u.uo 


1 7Q 4- n no 
1. lo ± u.uz 


l./o ± U.Uo 


1 7/1 4- n nQ 
1. / 4 it U.Uo 


i 7^ 4- n nQ 
l./o ± U.Uo 


1 qo _i_ n i /i 

l.oZ It U.14 






ft n*j 4_ n Qft 
O.Uo ± U.oo 


a n 4_ n o/i 
o.4y it U.z4 


O.zl ± U.zz 


^ Qn 4_ n q 
O.oU it U.oo 


c: 1 o 4- n oft 
O.lo ± U.zO 


c; /f q _L n p; i 
O.4o ± U.01 






i 7/1 4- n nft 
1. /4 ± U.UO 


1. (O it U.U4 


i 7/i J_ n nQ 

l.(4± U.Uo 


1 71 4_ n n/i 
1. / 1 it U.U4 


i 7/1 4- n nQ 
1. /4 ± U.Uo 


1 ot; j.n 11 
l.oO ± U. 11 






5 67 + 63 


5 33 + 24 


5 04 J- n 97 


5 05 + 25 


5 00 + 21 


5 47 _i_ n Qft 


F5blue 




1.72 ±0.10 


1.71 ±0.05 


1.68 ±0.04 


1.70 ± 0.03 


1.73 ±0.02 


1.79 ±0.06 






5.44 ±0.72 


5.06 ±0.27 


4.85 ± 0.24 


4.87 ±0.22 


4.79 ±0.18 


5.00 ±0.14 


F6blue 






1.75 ±0.20 


1.74 ±0.12 


1.75 ± 0.09 


1.83 ±0.06 


1.86 ±0.08 








5.27 ± 1.35 


4.90 ±0.44 


4.68 ±0.24 


4.71 ± 0.23 


4.83 ±0.14 


a = a = 


1.76 


^0 


ro 


ro 


ro 




ro 


Flbluc 




6.74 ±0.53 


6.02 ±0.15 


5.92 ± 0.18 


6.00 ±0.15 


6.00 ±0.21 




F2 b i uc 




6.23 ±0.14 


5.69 ±0.11 


5.44 ±0.19 


5.52 ±0.16 


5.39 ±0.15 


6.17 ±0.58 


F3bi uc 




5.98 ±0.12 


5.52 ±0.09 


5.25 ± 0.12 


5.32 ±0.15 


5.20 ±0.15 


5.34 ±0.14 


F4 b i uo 




5.70 ±0.16 


5.32 ±0.22 


5.06 ±0.18 


5.10 ±0.17 


5.03 ±0.18 


5.33 ±0.15 


F5bi uc 




5.51 ± 0.11 


5.14 ±0.11 


4.94 ±0.17 


4.93 ±0.13 


4.82 ± 0.09 


4.97 ±0.16 


F6bi uc 






5.29 ±0.73 


4.92 ± 0.27 


4.69 ±0.35 


4.62 ± 0.19 


4.68 ± 0.17 










TABLE 12 








The correlation length and slope between red and blue galaxies on small scales (0.16 h~ 


_1 Mpc < r p < 0.98 h 


_1 Mpc) 


Sample 




Fired 


F2 rc d 


F3 rcc i 


F4 rcd 


F5 re( j 


F6 rc d 






a 


Q 


a 


a 


a 


a 






ro 


ro 


ro 


ro 


ro 


ro 


Flbluc 




1.87 ±0.16 


2.01 ±0.13 


2.06 ± 0.09 


1.84 ±0.14 


i ns -U n 1 ^ 
jL.yo it u. io 








7.60 ±0.96 


5.08 ± 0.43 


4.90 ± 0.41 


6.14 ±0.56 


U . ( '4: ZL U.OO 




F 2 bluo 




2.04 ±0.17 


1.81 ±0.13 


1.74 ±0.08 


1.74 ±0.10 


1.0 1 it U. ID 


1 4- n i ft 
i.oy it u.id 






5.15 ±0.35 


5.31 ±0.35 


5.34 ±0.44 


5.50 ±0.36 


Qn i n a a 
O.oy ± U.44 


ft Cft 4- 1 Q(l 

o.oD it i.oy 


F3blue 




1.77 ±0.16 


1.77 ±0.12 


1.67 ±0.15 


1.67 ±0.07 


i k _l n nfi 
1.O0 ± U.UD 


i ftn 4_ n i o 
i.oy it U.lz 






6.03 ±0.53 


5.19 ±0.44 


5.27 ±0.34 


5.23 ± 0.48 




ft nn 4_ n c^n 
D.UU it u.oy 


F4 b l ue 




1.92 ±0.19 


1.68 ±0.15 


1.67 ±0.14 


1.67 ±0.11 


i ftyi _L n nc 
1.D4 ± U.Uo 


1 71 4_ n 1 o 
1. / 1 it U. IZ 






5.12 ±0.45 


5.36 ±0.56 


5.11 ± 0.68 


5.11 ±0.46 


5.04 ± 0.37 


5.82 ± 0.45 


F5blue 




1.85 ± 0.12 


1.84 ±0.16 


1.73 ± 0.08 


1.65 ±0.13 


1.56 ±0.06 


1.67 ±0.11 






5.28 ± 0.71 


4.71 ±0.58 


4.71 ±0.42 


4.89 ± 0.45 


5.09 ± 0.46 


5.16 ±0.48 


F6blue 






2.07 ±0.27 


1.79 ±0.24 


1.73 ±0.23 


1.68 ±0.11 


1.70 ±0.12 








4.29 ± 1.47 


4.92 ± 0.48 


4.66 ±0.27 


4.77 ±0.51 


5.04 ± 0.49 


a = a = 


1.76 


ro 


ro 


ro 


ro 


r 


ro 


Flbluc 




7.22 ± 0.61 


5.94 ±0.28 


5.98 ± 0.31 


6.53 ±0.36 


6.77 ±0.32 




F2bi uc 




6.26 ±0.31 


5.52 ±0.16 


5.25 ± 0.13 


5.43 ±0.16 


5.20 ±0.21 


6.02 ±0.55 


F3biu C 




6.06 ± 0.37 


5.22 ±0.13 


4.93 ± 0.12 


5.02 ±0.15 


4.97 ±0.15 


5.74 ±0.26 


F4 b i uc 




5.69 ±0.41 


5.11 ±0.35 


4.80 ± 0.34 


4.79 ±0.38 


4.75 ± 0.35 


5.61 ±0.28 


F5biu C 




5.59 ±0.27 


4.98 ±0.33 


4.61 ±0.25 


4.53 ±0.34 


4.44 ±0.19 


4.91 ±0.27 


F6biu C 






5.20 ±0.75 


5.01 ± 0.24 


4.57 ±0.19 


4.58 ±0.28 


4.90 ±0.24 



for red galaxies, and h/b* = 0.560 + O.IOIL^L* for blue 
galaxies. Note that the constant component in the rela- 
tion is much larger for red galaxies than for blue galaxies, 
although the linear terms are now similar. These results 
are in good agreement with those obtained in Zehavi et 
al. (2005) based on auto-correlation functions. 

6. NON-LINEAR AND STOCHASTIC BIAS 

In the last two sections we have analyzed the auto- 
and cross-correlation functions for galaxies of different 
luminosities and colors. We have seen that on scales 
r p > 1 /i _1 Mpc all correlation functions are roughly par- 
allel, but on smaller scales different galaxies may have 
different behaviors (see Fig. [5]). In order to investigate 
the scale-dependence of the correlation functions in more 



detail, we plot in Fig. [15] the ratios between the various 
projected correlation functions, W(r p ), and a normal- 
ization function, Wn(r p ), which is the power-law fit to 
the projected correlation function of the sample F* on 
large scales and corresponds to 7-n = 5.49 /i _1 Mpc and 
a = 1.76. Three top panels show results for red galax- 
ies; six middle panels show results between red and blue 
galaxies; and three bottom panels show results for blue 
galaxies. Based on the figure, we can draw the following 
conclusions: 

• On scales larger than 1 /i _1 Mpc, all the correla- 
tion functions are roughly parallel, a fact we have 
mentioned earlier. 

• The cross-correlation functions between the bright- 
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TABLE 13 

The average value of IZij over the large scales (0.98 fe _1 Mpc < r p < 9.6 7i _1 Mpc) 



Sample 


Fl 


F2 


F3 


F4 


F5 


F6 


Fl 




1.00 ±0.04 


1.05 ±0.06 


1.05 ±0.11 


1.14 ±0.08 




F2 


1.00 ± 0.04 




1.13 ± 0.06 


1.02 ±0.03 


1.12 ± 0.02 


0.91 ± 0.14 


F3 


1.05 ± 0.06 


1.13 ±0.06 




1.02 ±0.02 


1.20 ±0.06 


0.93 ±0.08 


F4 


1.05 ±0.11 


1.02 ±0.03 


1.02 ±0.02 




1.30 ±0.13 


0.93 ± 0.05 


F5 


1.14 ±0.08 


1.12 ±0.02 


1.20 ±0.06 


1.30 ±0.13 




0.94 ± 0.02 


F6 




0.91 ±0.14 


0.93 ±0.08 


0.93 ±0.05 


0.94 ±0.02 






Fl rod 






F4 rcd 


F5 re d 


F6 rcd 


Fl rcd 




0.98 ±0.05 


0.98 ±0.05 


1.05 ±0.19 


1.09 ± 0.12 




F2 rcd 


0.98 ± 0.05 




1.02 ± 0.05 


1.04 ± 0.06 


1.06 ± 0.05 


1.17 ±0.27 


F3 rcd 


0.98 ± 0.05 


1.02 ±0.05 




0.99 ± 0.03 


0.99 ± 0.07 


0.98 ±0.13 


F4 rcd 


1.05 ±0.19 


1.04 ±0.06 


0.99 ±0.03 




0.97 ±0.04 


1.03 ± 0.08 


F5 rcd 


1.09 ±0.12 


1.06 ±0.05 


0.99 ± 0.07 


0.97 ±0.04 




0.94 ± 0.03 


F6 rcd 




1.17 ±0.27 


0.98 ±0.13 


1.03 ±0.08 


0.94 ± 0.03 






Flbluc 


F2 b iuc 


F3 b iuc 


F4 b iuc 


F5 b iuc 


F6 b iuc 


Flbluc 




1.20 ±0.23 


1.41 ±0.28 


1.76 ±0.51 


1.94 ±0.71 




F2 b i uc 


1.20 ±0.23 




1.06 ±0.09 


1.10 ±0.19 


1.14 ± 0.21 


1.45 ±0.31 


F3biuc 


1.41 ±0.28 


1.06 ±0.09 




1.01 ±0.07 


1.08 ±0.10 


1.10 ±0.10 


F4 bluc 


1.76 ±0.51 


1.10 ±0.19 


1.01 ±0.07 




0.99 ±0.05 


1.01 ±0.05 


F5 b l ue 


1.94 ±0.71 


1.14 ±0.21 


1.08 ±0.10 


0.99 ± 0.05 




0.98 ± 0.04 


F6bi uc 




1.45 ±0.31 


1.10 ±0.10 


1.01 ±0.05 


0.98 ± 0.04 






Flbluc 


F2biuc 


F3biuc 


F4 b iuc 


F5biuc 


F6biuc 


Fl red 


1.31 ±0.12 


1.22 ±0.20 


1.10 ±0.13 


1.19 ±0.16 


1.14 ±0.17 




F2 rcd 


1.49 ± 0.40 


1.18 ±0.21 


1.10 ±0.11 


1.08 ±0.08 


1.11 ± 0.06 


1.00 ±0.15 


F3 red 


1.48 ±0.50 


1.21 ±0.24 


1.18 ±0.15 


1.12 ±0.08 


1.13 ± 0.21 


0.98 ±0.14 


F4 rcd 


2.12 ± 1.10 


1.35 ±0.28 


1.18 ±0.17 


1.16 ±0.15 


1.10 ±0.15 


1.00 ±0.16 


F5 rcd 


2.05 ±0.65 


1.42 ±0.39 


1.31 ±0.29 


1.18 ±0.20 


1.15 ± 0.19 


1.06 ±0.17 


F6 red 




2.27 ±0.96 


1.63 ±0.50 


1.33 ±0.27 


1.33 ± 0.40 


1.31 ±0.28 



TABLE 14 

The average value of IZij over the small scales (0.16 fe _1 Mpc <r p < 0.98 /i _1 Mpc) 



Sample 


Fl 


F2 


F3 


F4 


F5 


F6 


Fl 




0.92 ±0.10 


0.92 ± 0.08 


0.80 ± 0.34 


0.81 ± 0.44 




F2 


0.92 ±0.10 




1.07 ±0.11 


1.00 ±0.07 


1.23 ±0.23 


1.03 ±0.53 


F3 


0.92 ± 0.08 


1.07 ±0.11 




1.05 ±0.06 


1.13 ± 0.12 


0.84 ±0.06 


F4 


0.80 ± 0.34 


1.00 ±0.07 


1.05 ±0.06 




1.12 ±0.07 


0.90 ±0.11 


F5 


0.81 ±0.44 


1.23 ±0.23 


1.13 ±0.12 


1.12 ±0.07 




0.95 ± 0.04 


F6 




1.03 ±0.53 


0.84 ± 0.06 


0.90 ±0.11 


0.95 ± 0.04 






Flrcd 


F2 rcd 


F3 ro d 


F4 rcd 


F5 ro d 


F6rcd 


Flrcd 




0.84 ±0.11 


0.78 ±0.16 


0.73 ±0.20 


1.09 ± 0.38 




F2 rc d 


0.84 ±0.11 




1.03 ±0.13 


1.06 ±0.08 


1.20 ±0.15 


1.29 ±0.81 


F3 rc d 


0.78 ±0.16 


1.03 ±0.13 




1.06 ±0.09 


1.08 ± 0.14 


1.06 ± 0.32 


F4 rcd 


0.73 ±0.20 


1.06 ±0.08 


1.06 ±0.09 




1.06 ±0.09 


1.02 ± 0.08 


F5 rc d 


1.09 ±0.38 


1.20 ±0.15 


1.08 ±0.14 


1.06 ±0.09 




1.01 ±0.08 


F6 rcd 




1.29 ±0.81 


1.06 ±0.32 


1.02 ±0.08 


1.01 ± 0.08 






Flbluc 


F2 b iuc 


F3biuc 


F4 b iuc 


F5 b iuc 


F6biuc 


Flbluc 




0.60 ±0.62 


0.95 ±0.76 


0.79 ±1.08 


1.05 ± 1.58 




F2 b iuc 


0.60 ±0.62 




1.02 ±0.22 


1.10 ±0.19 


1.11 ±0.56 


0.63 ±0.73 


F3biuc 


0.95 ± 0.76 


1.02 ±0.22 




1.11 ±0.18 


1.21 ± 0.18 


1.13 ±0.30 


F4 b i uc 


0.79 ± 1.08 


1.10 ±0.19 


1.11 ±0.18 




1.00 ±0.11 


0.93 ± 0.22 


F5blue 


1.05 ± 1.58 


1.11 ±0.56 


1.21 ±0.18 


1.00 ±0.11 




1.00 ±0.09 


F6biuc 




0.63 ±0.73 


1.13 ±0.30 


0.93 ±0.22 


1.00 ±0.09 






Flbluc 


F2biuc 


F3biuc 


F4 b iuc 


F5biuc 


F6 b iuc 


Flrcd 


1.12 ±0.61 


1.58 ±0.60 


1.59 ±0.57 


1.54 ±0.35 


1.43 ± 0.59 




F2 red 


1.07 ±0.99 


1.41 ±0.25 


1.33 ±0.19 


1.28 ±0.26 


1.26 ±0.17 


0.94 ±0.60 


F3 re d 


0.94 ± 1.15 


1.63 ±0.37 


1.51 ±0.23 


1.43 ±0.19 


1.47 ±0.17 


1.12 ±0.08 


F4 rcd 


1.46 ± 1.94 


1.74 ±0.25 


1.75 ±0.27 


1.60 ±0.20 


1.57 ±0.16 


1.20 ±0.24 


F5 rcd 


1.28 ± 1.62 


2.91 ± 1.35 


2.42 ±0.57 


2.13 ±0.58 


2.11 ±0.52 


1.62 ±0.18 


F6rcd 




1.85 ±1.72 


3.40 ± 1.55 


2.68 ±0.68 


2.90 ±0.83 


2.43 ± 0.63 



The cross-correlation between galaxies 
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Fig. 12. — The relative bias parameters for red galaxies obtained 
from the correlation functions over 0.98 /i _1 Mpc to 9.8 /i _1 Mpc 
using 13 different approaches (see the main text for details). The 
error-bars are based on the scatter of the data points in the range 
0.98 /i _1 Mpc < r p < 9.8 h _1 Mpc. Panels from top to bottom 
correspond to the bias parameters for Fi re( j (i = 1,2,3,5,6), nor- 
malized by that of F4 ro[ j. The horizontal axis labels the different 
approaches that are used to determine the bias parameter for each 
subsamples. Some data points are missing, because the redshift 
overlap between Fl and F6 is too small to measure their cross- 
correlation function reliably. 

est red galaxies and other red galaxies have 
enhanced clustering on scales smaller than ~ 
1 /i _1 Mpc, and the enhancement is the strongest 
for the cross-correlation between the brightest and 
faintest samples. 

• The cross-correlation functions among faint blue 
galaxies, and between faint-blue and faint-red 
galaxies, appear to be suppressed slightly on scales 
< 1 /i _1 Mpc relative to the extrapolations from 
larger scales. 

• The brightest blue galaxies show enhanced corre- 
lation with red galaxies of different luminosities on 
scales < 1 /i _1 Mpc. 

Some of these results are expected, as it is well known 
that red galaxies reside preferentially in rich systems, 
such as clusters and groups, while faint blue galax- 
ies are distributed in the 'field' (i.e., in dark matter 
haloes of significantly lower mass). Note, however, that 
the cross-correlation functions between the brightest red 
galaxies and faint red galaxies are lower than the auto- 
correlation function of the brightest galaxies on scales 



Fig. 13. — The same as Fig. 1121 for blue galaxies. Panels from 
top to bottom correspond to the bias parameters for Fibi uc (i = 
1, 2, 3, 5, 6), normalized by that of F4bi uc . 

r P ^ 1 h~ l Mpc, suggesting that not all faint red galax- 
ies are associated with the brightest red galaxies. There 
is no enhancement in the cross-correlation functions on 
small scales between red galaxies and all except the 
brightest blue galaxies. This again can be explained by 
the fact that blue galaxies with low and intermediate lu- 
minosities are distributed in the field. 

One surprising result is that the brightest blue galax- 
ies are strongly correlated with red galaxies on small 
scales, and in many ways the clustering properties of 
these galaxies are similar to that of red galaxies with 
intermediate luminosities. In order to understand this 
result in more detail, we have examined the properties 
of this population of galaxies more closely. Using the 
group catalogue constructed by Yang et al. (2007, in 
preparation), we find that ~ 15% (~ 80%) of the bright- 
est blue galaxies are the central galaxies of groups with 
masses > 10 13 - 5 /i _1 M Q (> 10 13 /i _1 M Q ), compared to 
~ 30% (<~ 85%) for the red galaxies in the same lumi- 
nosity bin. It is therefore possible that a fraction of the 
brightest blue galaxies are in fact 'early-type' galaxies 
and their relatively blue colors are probably due to re- 
juvenation of star formation or AGN activities in their 
centers. Clearly, further analysis is required in order to 
understand the properties of this population of galaxies 
in more detail. 

Next we discuss how our results may put constraints on 
the stochasticity /non-linearity in the bias relation. Con- 
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logioVk log 10 Li/L* 

Fig. 14. — Left panels: The average relative bias parameter, obtained using the correlation data in 0.98 h~ 1 Mpc < r p < 9.8 h~ 1 M.pc 
as a function of galaxy luminosity, for red (upper panel) and blue (lower panel) galaxies. The results are obtained using the flux-limited 
samples, and both the bias parameters and the luminosities are normalized by the values ofF4 rc( j and F4t,i uc , respectively. For a given 
luminosity, the mean and error-bar are based on the different measurements shown in Figs. 1121 and |13I The solid curves show fits to the 
data. Right panels: The same as the left panels, except that here the bias parameters and the luminosities are normalized by the values 
for L* galaxies. 



sider the density field (smoothed on some scale) traced 
by a population i of galaxies, Si. In general we may write 
it in terms of the density field of another population, Sj , 

as 



Si — BijSj 



(13) 



3 J C «J ) 

where the first term on the right-hand side describes the 
linear part of the correlation between the two popula- 
tions, the second term is the non-linear part, and is 
the stochastic, uncorrelated part: (tijSj) = 0. It then 
follows that 



(S i S j )=B ij (S^ + (S j f i (S j )), 



(14) 



and 

{5?) = S?.(5]> + (4) + </?&)) + 2i?« <Wi)> • (15) 

Thus, for any two populations we can define two relative 
bias factors, one based on the cross correlation with an- 
other population, and one based on the auto-correlation 



functions of the two populations: 

? (c) _ (8j6j) 



D 



n - /x2\ ' 



B 



w = mi 



1/2 



(16) 



(17) 



The square of the ratio between these two bias factors is 




1 + 



(4-) + & 
<«7<*?> ' 



where 



(18) 



(19) 
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Fig. 15. — The ratios between the projected correlation function W(r p ) and Wo(r p ) for cases with various combinations of color and 
luminosity, where Wo(r p ) is the power-law fit to the projected correlation function of sample F* on large scales (ro = 5.49 ft _1 Mpc, 
a = 1.76). Three top panels show the results for red galaxies, six middle panels show the results between red and blue galaxies, and three 
bottom panels show the results for blue galaxies. In each panel, iR-jR denotes Fi re( j with Fj re( j, iB-jB denotes Fibi uc with Fjbi ue > iR-jB 
denotes Fi rct j with Fjt,i ue - 



describes the contribution due to the non-linear part of 
the relative bias relation. Note that TZ is related to the 
linear correlation coefficient, r, defined in Dekel & Lahav 
(1999; their eq.[13]) as TZ — r~ 2 . For two populations, i 
and j that are entirely uncorrelated one has that Rij — > 
oo, while Rij = 1 when i and j are perfectly correlated. 
Both the stochastic component, (e?-), and the non-linear 
component, £y, can cause TZij to deviate from unity, and 
in general it is difficult to separate stochasticity from 
non-linearity. On large scales, however, the non-linear 
component in the bias relation may be neglected, and we 
have that 



7? — 



(4) 



(20) 



To examine the importance of stochasticity /non- 
linearity in the relative bias relations, we calculate TZij 
for galaxies of different luminosities and colors. Based 
on equation (TTg|) . we estimate TZij using 



TZ; 



WgWj. 
W 2 



(21) 



where Wu and Wjj are the auto-correlation functions of 
subsamples 'j' and 'j', respectively, and Wij is the cross- 
correlation function between the two subsamples. Note 
that H' and 'j' denote not only luminosity samples but 
also color samples. Unlike what we did above, here we 
estimate Wu, Wjj and all in the overlapping region 
of the two samples in question, so that volume effects are 
minimized. 
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Fig. 16. — The value of TZij as a function of r p for cases with various combinations of color and luminosity. Panels from top to bottom 
show results for all-all, red-red, blue- blue, and red-blue samples, respectively. The straight line in each panel shows IZjj = 1. 



In Fig llOl wo show IZjj as a function of r p for various 
combinations of color- and luminosity-samples. The re- 
sults are quite noisy in some cases, again because of the 
small overlap of the samples involved. As one can see, for 
cases where TZij is well determined, it is close to 1 for all- 
all and red-red pairs of samples. For blue-blue pairs, TZ^ 
is also close to 1 except the cases involving the bright- 
est galaxies where at large scales TZij can be significantly 
larger than 1. As we have discussed above, the bright- 
est blue galaxies have clustering properties different from 
that of fainter blue galaxies, because a significant frac- 
tion of them are, like bright red galaxies, central galaxies 
of relatively massive groups. For red versus blue galax- 
ies, TZij is significantly larger than 1 on small scales. The 
signal is also more prominent between faint red and faint 
blue galaxies and can extend to a few /i _1 Mpc. 

To show this more clearly, we plot the average value 
of TZij over two different scales. The results are shown 
in Figs. OH -[Hand listed in Tables [13] -[H for all ver- 
sus all, red versus red, blue versus blue, and red versus 



blue, respectively. In all the figures, circles show results 
based on clustering on large scales (0.98 /i -1 Mpc < r p < 
9.6/i _1 Mpc), while crosses show results based on cluster- 
ing on small scales (0.16 /i _1 Mpc < r p < 0.98 /i _1 Mpc). 
As one can see, on large scales (> 1 h~ 1 Mpc) TZij is very 
close to 1 for red-red, blue-blue (with the exception of 
cases involving the brightest blue galaxies), and all-all 
pairs, but has a value between 1.4 to 2.2 for pairs be- 
tween faint red and blue galaxies. On scales < l/i _1 Mpc, 
the ratio is again very close to 1 for all-all, red-red and 
blue-blue pairs, but is significantly larger than 1 for red 
versus blue galaxies and is as large as ~ 3 for faint-red 
versus faint-blue galaxies. 

These results suggest that red (or blue) galaxies of 
different luminosities trace each other almost linearly 
and in an almost deterministic way on both small and 
large scales, but there is a significant stochastic/non- 
linear component in the relationship between red and 
blue galaxies of different luminosities, and in particular 
on small scales. The large value of Rij on small scales 
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Fig. 17. — The ratio IZij defined in equation ( 1181 for all (left panels) and red (right panels) gala 



In all panels, circles show results 



based on data on large scales, 0.98 /i 1 Mpc < r p < 9.6?t 1 Mpc, while crosses show results based data on small scales, 0.16/i 1 Mpc < r p < 
0.98 ?t — 1 Mpc. From top to bottom are results for i = 1,2,3,4,5,6, respectively. The value of j for each case is labeled on the horizontal 
axis. The dashed line in each panel shows IZij = 1. The error-bars are estimated using the scatter of the data points on large scales 
0.98 h _1 Mpc < r p < 9.6 /i _1 Mpc and small scales 0.16 h~ 1 Mpc < r p < 0.98 /i _1 Mpc, respectively. 



between red and blue galaxies (except the brightest blue 
galaxies) is almost certainly due to the fact that blue 
galaxies of low and intermediate luminosities avoid dense 
groups and clusters while red galaxies reside preferen- 
tially in such places. The fact that such stochastic/non- 
linear component extends to a few /i _1 Mpc (see Fig. ITB)) 
suggests that the clustering of clusters and groups on 
large scales may produce segregation of the galaxy pop- 
ulation on large scales. 

7. SUMMARY 

In this paper, we use galaxy samples constructed from 
the SDSS Data Release 4 (DR4) to study the cross- 
correlation functions between galaxies with different lu- 
minosities and colors. Galaxies are separated at g — r = 
—0.7 — 0.032(^/^0.1 + 16.5) into red and blue samples 
according to their colors, and each of the color samples 
is further divided into 6 subsamples according to lumi- 
nosity. We measure the projected two-point correlation 
function for each subsample, and the projected cross- 
correlation for each pair of subsamples. Each correla- 



tion function is fitted with power laws over two different 
ranges of separations: 0.98 ft^Mpc < r p < 9.6 /i _1 Mpc 
and 0.16 fi^Mpc < r p < 0.98 /i _1 Mpc, and each of the 
power laws are characterized by a correlation length, ro, 
and a logarithmic slope, a. Our main results can be 
summarized as follows: 

1. At projected separations r p > 1 /i _1 Mpc, all corre- 
lation functions are roughly parallel to each other 
(but see exceptions for the brightest galaxies with 
-22.0 > M r ,o.i > -23.0 in Zehavi et al. 2005; Li 
ct al. 2006), but on smaller scales different popu- 
lations of galaxies have different behaviors. 

2. The auto- and cross-correlation functions of red 
galaxies on r p < 1 /i _1 Mpc are significantly en- 
hanced relative to those on larger scales. The effect 
is the strongest between the brightest red galaxies 
and other red galaxies. Such enhancement is ab- 
sent for blue galaxies and in the cross-correlation 
between red and blue galaxies. 
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Fig. 18. — The same as Fig. 1171 but here results are shown for blue versus blue galaxies (left panels) and red versus blue galaxies (right 
panels). 



3. For blue galaxies the luminosity-dependence of the 
correlation amplitude on large scales is strong over 
the entire luminosity range, while for red galaxies 
the dependence is weaker and becomes insignificant 
for luminosities below L* . There is an indication 
that the luminosity-dependence may inverse at the 
faint end. 

4. For a given luminosity, red galaxies are more 
strongly clustered than blue galaxies on both large 
and small scales, and the difference is larger for 
fainter luminosities. 

5. One large scale, the bias factors of a given popula- 
tion of galaxies obtained from its cross correlations 
with other populations of galaxies are similar, and 
comparable to that based on the auto-correlation 
function. 

6. On both large and small scales, the ratio TZij is very 
close to 1 for red (or blue except the brightest) 
galaxies of different luminosities, suggesting that 
the bias relations between red (or blue) galaxies of 
different luminosities are close to linear and deter- 
ministic. 



7. On scales < 1 /i _1 Mpc, Rij is significantly larger 
than 1 between red and blue galaxies and as large 
as ~ 2 between faint red and faint blue galaxies, 
suggesting that a significant stochastic/non-linear 
component exists in the bias relations between blue 
and red galaxies on small scales. 

8. The clustering properties of the brightest blue 
galaxies resemble more that of red galaxies of in- 
termediate luminosities than that of fainter blue 
galaxies, indicating that a significant fraction of 
them are located in relative rich systems. 

The results obtained in this paper provide important 
information about the relationships between the spatial 
distributions of galaxies of different properties in space. 
Part of our findings are consistent with previous probes 
(e.g. the luminosity and color dependences of SDSS 
galaxies based on the auto correlation functions by Ze- 
havi et al. 2005; Li et al. 2006). Most of our findings will 
put new constraints on the galaxy formation models. To 
explore the implications of our findings for galaxy forma- 
tion in the cosmic density field, the results obtained here 
should be combined with semi- analytical and/or halo oc- 
cupation models to constrain how galaxies populate dark 
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matter halos, and how galaxies are distributed in indi- 
vidual dark matter halos. We will come back to this in 
the future. 

During the final stages of this project a paper appeared 
by Swanson et al. (2007) , who carried out a similar inves- 
tigation to study the relative biases for galaxies of dif- 
ferent luminosities and colors. Although their analysis 
is based on counts-in-cells, as opposed to the correlation 
functions used here, their results are qualitatively similar 
to ours. A detailed comparison, however, is complicated 
because of the differences in the statistical methods and 
sample selections. 
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